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ABSTRACT 

In this paper we consider the laminar boundary layer which forms on a sharp cone in a 
supersonic freestream, where lateral curvature plays a key role in the physics of the problem. 

This flow is then analysed from the point of view of linear, temporal, inviscid stability. 
Indeed, the basic, non- axi symmetric disturbance equations are derived for general flows of 
this class, and a so called “triply generalised” inflexion condition is found for the existence of 
“subsonic” neutral modes of instability. This condition is analogous to the well-known gen- 
eralised inflexion condition found in planar flows, although in the present case the condition 
depends on both axial and aximuthal wavenumbers. 

Extensive numerical results are presented for the stability problem at a freestream Mach 
number of 3.8, for a range of streamwise locations. These results reveal that a new mode of 
instability may occur, peculiar to flows of this type involving lateral curvature. 

Additionally, asymptotic analyses valid close to the tip of the cone / far downstream of 
the cone are presented, and these give a partial (asymptotic) description of this additional 
mode of instability. 


1 This research was supported by the National Aeronautics and Space Administration under NASA Con- 
tract No. NAS1-18605 while the second author was in residence at the Institute for Computer Applications 
in Science and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 23665. 
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.L Introduction and motivation 

It is well known that unlike many of their incompressible counterparts, 
supersonic boundary layers arc generally susceptible to inviscid forms of 
instability; this comment, in particular, applies to Blasius-lype boundary 
layers. The paper which first tackled the problem of stability of 
supersonic boundary layers, from any kind of rigorous mathematical 
standpoint was that of Lees and Lin (1946), who presented some extremely 
important results on the problem, including the importance of the quantity 


■v ♦ 

— „ [p* — „] (where u* denotes velocity tangential to the surface, y* 
dy 1 dy J 

the normal to the surface, and p* the fluid density), which plays the 
same role as that of an inflexion point in incompressible flows. In 
particular at the point where the above expression is zero (y* = y 0 *, say), 
then a neutral mode exists with wavespeed u*(y 0 *) (although see the 
restriction (4.21) later in this paper). Lees and Lin also showed how 
viscosity played a similar role at a critical layer as in analogous 
incompressible flow situations. This two-dimensional work was later 
extended to three-dimensions by Reshotko (1962). 

The development of digital computers enabled accurate numerical 
computations of the problem to be made. One of the earliest was that of 
Brown (1962), and this was followed by a number of studies by Mack 
(1963, 1964, 1965a, b, 1969, 1984, 1987a), both of an inviscid (large 
Reynolds number) and viscous (finite Reynolds number) nature. These 
computations revealed a further distinction between supersonic and 
incompressible flows, namely that in the latter case an infinite sequence of 
modes is possible. These are important, because generally it is one of 
these "higher" modes that exhibits the largest growth rate according to 
inviscid, linear stability theory. 

In spite of their great practical importance, flows of this type, but 
involving lateral curvature have received scanty attention over the years. 
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Duck and Hull (1989a) showed how in supersonic flows, curvature interacting 
with viscosity could provoke additional instabilities (axi symmet ric in 
form), provided the body radius was below some critical value. Duck and 
Hall (1989b) then went on to show how a similar effect occured with 
non- axi symmet ri c modes (which, in fact, turn out to be generally more 
unstable than corresponding axisymmetric modes). 

Recently, Duck (1989) (hereafter referred to as I) has shown the 
influence of curvature on the axisymmetric inviscid stability of supersonic 
boundary layers, in particular those that form on thin straight circular 
cylinders. It was found that curvature has a profoundly stabilising effect 
on these modes. Interestingly, it was shown how curvature altered the 
generalised inflexion condition described above, and a modified (or "doubly 
generalised") inflexion condition, involving the radius of curvature was 
derived . 

In the present paper we extend the work of I to non-axi symmet ric- 
disturbances which, indeed, turn out to be more important than axisymmetric 
disturbances considered previously. Further, rather than studying/applying 
our techniques to the thin straight circular cylinder, we consider a 
somewhat more practical configuration, namely that of a sharp cope. 

To (considerably) facilitate already lengthy and complex computations, 
following I we assume that the tip of the cone has a finite radius, and 
that the associated boundary layer on the cone surface is planar at the cone 
tip. (We expect that far downstream the flow will become planar once again, 
having utilised the Mangier, 1948, transformation.) Following previous 
studies in this area, we also ignore the presence of any shocks. 

Mack (1987b) has performed some computations for the stability of the 
flow over a cone in supersonic flow, at finite Reynolds numbers, but found 
little difference with corresponding planar results. Here, we deliberately 
allow curvature to occur throughout the study, both in the equations 
governing the basic flow, and in the disturbance equations. 
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A few preliminary numerical results of this study have been presented by 
Duck (1990). Here, we present additional numerical results, and also a 
number of important asymptotic results related to these calculations. 


2. Equations of motion and state 

The general layout of the problem is shown in Fig. 1. The z* axis 
lies along the cone axis, r* denotes the radial coordinate, and 0 the 
azimuthal coordinate. The velocity vector v* has components vj*, V2 
and vj* in the r*. 0 and z* directions respectively. Although we 
shall be concerned with a basic flow which is independent of 0 (and with 
v 2 * = 0), when we go on to consider the stability of the flow, we shall be 
concerned with non-axisymmetric disturbances. 

In the cylindrical polar coordinate system as defined above, the full 
equations of continuity, momentum and energy take on the following forms 


(Thompson 1972) 


K * X CP* v,*) ♦ i. L (P* v 2 *> 
dt or r o0 


° f * * \ 

+ ^*(P v 3 ) 


* * 

+ ? 4 i = o. 


>*[ 


Dvj* (v 2 *) 2 
Dt* ’ "T* 


dp 3 S rV 
~ ' 37* + 3r 


+ iy + d h* zL 

r 30 N dz 

+ ^r*r* ~ ^66 f 
r 

* r dv/ 1 = . ' §e*+ 1 

p lor J 7 * 30 Hr 


+ i_ 3 -96 * 3 Uz * 
7* 30 ' 3z* 


( 2 . 1 ) 


( 2 . 2 ) 


+ 2 lr*0 

i • 


(2.3) 
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p * *£ , . §e! + LLCil + > a 

Dt 3z dr r 30 


+ - -|v ^ + . 

dz r 


(2.4) 


p ' 57 <c p 


t*) . 2e^ = r* + I,, i-y [ kV 1 

Dt 7 r 7 dr* l 3r* J 

+ * L ( K* 2C | 

( r) 2 30 l 30 ' 

3 r K*3T* | 

+ ^ l 3^ I' 


(2.5) 


Here p* is the density of the fluid, p* the pressure, Cp the 
specific heat (at constant pressure), and K* the coefficient of heat 
conduction. The Eulcrian operator is defined as 


D = 3_ 
Dt T 3t 


_ _ vo* 3 * 

_ + Vl + -f- - + v 3 


* L 

3r* r* 30 


3 

3? 


(2.5) 


and the viscous stress components, assuming Newtonian flow are defined 
to be 


x r V = V -^4 + x* v, v* , 

dr 

Xee - V [ ft + 7 ^ ] + >■* v I* . 


IjV ■ 2H* + X* 7.v* . 

a z 


I,*. - l.,W [ i. L, [: 

V f. W . ?[!$.!&]. 

The dispersion function f* in (2.5) is given by 



( 2 . 6 ) 

(2.7) ' 

( 2 . 8 ) 

(2.9) 

( 2 . 10 ) 
(2.11) 
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r* = 2ji |d 2 v + d§ 0 + d 2 Y 


+ 2 D?*e + 2 D^0 + 2 D^V 


+ (X* - 2 JI*) (V.V*) 2 . 

I 


( 2 . 12 ) 


Here the 'D' terms are the components of the rate-of-deformation tensor, 
for example 


n 1 ■ vi* 

D ee — *■ —£r + ~r~ ■ 


n * 1 

D Z 0 = * 


r 30 r 
* 


3v 2 * i 3v_ 

dz + 1 * ae 


(2.13) 


(2.14) 


The coefficients p* and X* above denote the first coefficient of 
viscosity and bulk viscosity respectively (which are expected to be 
functions of temperature). 

The equation of state is taken to be that pertaining to a perfect gas, 

i .e. 


p* = p* R* T*. (2.15) 

where R* is the gas constant. 

With reference to Fig. 1, the surface of the cone is taken to lie 
along r* = a* + Xi z* , z* > 0, (later, important assumptions regarding 
the size of the slope parameter will be made), and so on this surface we 
requi re 

vj* = v 2 * = v 3 * = 0. (2.16) 

If the surface of the cone is insulated then the following additional 
boundary condition must be imposed 

= 0 . (2.17) 

(where n* denotes an outwards normal to the wall). In the case of heated/ 
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cooled walls, then the condition 


* * 

T = T 

w 


(2.18) 


must be imposed at the surface. 

Conditions remain to be specified at z* = 0; for this we follow I 
(precisely) by assuming that the boundary layer at this location has zero 
thickness, enabling planar conditions to be imposed at this position. 
Assuming the cone to be slender, then the far-field conditions are taken to 


be uni form, with 


* * „ 

vj = V2 = 0 

v 3 * = U„* 


T* = T * 


(2.19) 

( 2 . 20 ) 
( 2 . 21 ) 


(Indeed, in the case of non* slender cones, if the shock wave is attached, 
downstream of the shock the external flow is i rrotat ional , and there 
is also a constant slip velocity at the surface of the cone.) We next go on 
to derive the basic (boundary- 1 ayer) flow on the surface of the cone, 
assuming curvature plays a key role in the physics of the problem. 


We define our Reynolds number on the tip radius of the cone, a , as 


follows 


Re = LL* a* p w * 


(3.1) 


and this will be taken to be large throughout this paper. As noted 
previously, the basic flow is taken to be independent of 6 , and has no 
azimuthal velocity component (i.e. V 2 * = . 0 ) . 

It is now convenient to introduce non-dimensional parameters as follows 
(v j , v 3 , r, z, T, p, p) = 


f „ * * 

[Re V] V 3 


* * 
r z 


^ t — 7 t 3 

Ugo a Rea 


(3.2) 


ii if 1 1 i'll ii 
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A key element of this paper (as in I) is the inclusion of curvature 


terms in the governing equations to leading order. If this is'lo be the 
ease, then we must have that 

X] = Re" * X 

where X = 0(1) , (3.3) 

implying a slender cone. 

The leading order governing equations may then be written (assuming 
Re — > «) 

57 PH + ?i * H R) = 0 <3 ' 4) 

3 ?- 0 < 3 - 5 ) 

3 ? *«3 p ■ ? It h If 1 ] < 3 « 

V 137 * v 3 g - MT(r-l) M. 2 [|^] 2 


t a [rp ar 

+ 7 IT [5T JF 


where the result 


P = 


has been used, o is the Prandtl number, namely 


M- c p 

o = J- 

K 


(3.7) 


(3.8) 


(3.9) 


(which is taken to be a constant in this paper), y denoting the ratio of 
specific heats, and is the freest ream Mach number, namely 


M«, = Uo»* / (y R* Tj)i . 


(3.10) 


The boundary conditions to be applied to this system are 
v 1 = v 3 = 0 on r = 1 + Xz , 

V3 — * 1 , T — » 1 as r -» » , (3.11) 

together with a wall temperature condition; in the. case of insulated walls, 
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(lo leading order) 

= o on r = 1 + Xz , (3.12) 

whilst in the case of heatcd/coolcd surfaces 

T = T w on r = 1 + Xz . (3.13) 

To close the problem, a viscosi ty/ temperature law is required. For the 
purposes of this paper we take the linear Chapman law (Stewartson 1964), 
namely 

ji. = CT (3.14) 

where C is taken to be constant (although as noted in I, there would be 
no conceptual difficulties in taking more complex variations of viscosity 
with temperature). 

If we write 
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As described previously we assume planar conditions prevail at z = 0, 
where the boundary layer is taken to have zero thickness. The problem is 
thus singular at z = 0, and consequently sealed variables must be 
introduced in order to (numerically) solve the system (3.16) - (3.19) 
accurately. Specifically we write 


vi = C ' 1 vj (n.O. 

V3 ='v 3 Ot.O. 

T-Tdi.O. (3.20) 

where 

C = 2* . (3.21) 

and Tj = (r - 1 - (3.22) 

The quantities with a hat are then expected to behave regularly as 
£ — » 0, and indeed approach the planar solution. Equations (3.16) - 
(3.19) now take the following form 


d 





(3.23) 



(3.24) 




£v3 d\'. 


v 3 



(3.25) 


* 9T Cv 3 9T 
V 1 3Ti * ■ 


v 3 



A 

9T 


= T2(y-1)M„2 


f - 1 2 
9v 3 T9 

frT 9T1 

9q + r 

[o~ 3iJ 


(3.26) 
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Here of course 

r = 1 + \(, 2 + CO • (3.27) 

The boundary conditions in terms of these variables arc 
v l = v 3 = 0 on T| = 0 

V 3 — ♦ 1 , T — * 1 as !]—»“>. (3.28) 

In the case of insulated walls, the additional surface condition is 

|X = 0 on v\ = 0 . (3.29) 

whilst for heated/cooled walls, 

T = T w on r| = 0 . (3.30) 

Setting £ = 0 reduces the system to an ordinary differential system 
(corresponding to the planar case) as in I in the same limit (indeed, 
setting X = 0 reduces (3.23) - (3.26) to the corresponding system 
considered in I). The solution of the ordinary differential system at 
£ = 0 then provides the initial conditions for a (straightforward) 
Crank-Nicolson scheme in identical to that used in I. For the 

purposes of this paper, we shall focus our attention on the insulated wall 
case (although the heated/cooled case may be treated in exactly the same 
manner). Distributions of wall temperature with axial coordinate £(= fi) 
are shown in Fig. 2a (M* = 2.8) and Fig. 2b (M„ = 3.8), and the 
corresponding distributions of wall shear v^Iti = 0 are shown in Fig. 3a 
(M*. = 2.8) and Fig. 3b (M* = 3.8). 

In all cases, these distributions are quite different to the 
corresponding X = 0 distributions (i.e. the distributions on a straight 
circular cylinder) as found in I, where at both Sfc = 2.8 and 3.8 
a monotonic decrease in values was found. It is also quite clear that the 
results evolve from the planar case, to the far downstream limit, as 
predicted by the Mangier Transformation (Mangier 1946, Stewartson 1964), 
namely 


11 


V 3 t i|ti = o. ; 
T U = o, c - 


^ 3 ” v 3n It, = 0. c = 0 


(3.31) 


n = o, ; = o 


In the following section we go on to investigate the stability of 
flows of this class, subject to small amplitude inviscid disturbances. 


4. Inviscid stability of the flow 
4.1 Distu rbance equations 

In this section we derive the disturbance equations relevant to small 
amplitude disturbances in anv supersonic axi symmetric boundary layer type 
flow. In I, just axisymmetric disturbances were considered; here we 
consider the more general case of non-axi symmetric perturbations of the 
flow. 

We consider disturbances whose wavelength in the axial direction 
(1/a) is comparable to the (tip) radius of the cone. Specifically, at a 
fixed z station we write 

vj* = 5 ft v, (r) E + 0(5 2 ) . 

v 2 * = 5 U„* v 2 (r) E + 0(5 2 ) , 

v 3 * = U M * [W 0 (r) + 5 v 3 (r) E] + 0(8 2 ) , 

T* = [T 0 (r) + 5 T (r) E] + 0(5 2 ) , 

p = Pm + 8p (r) e] + 0(5 2 ) , 

p = p„* R* T„* [1 + 5p (r) E] + 0(8 2 ) , (4.1) 

where E = exp [i a (z - ct) + in 9J , (4.2) 

and 8 is the scale of the disturbance (taken to be dimini shingly small), 
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whilst 

i = (Oa*) t* . 
z = z*/a* 

W 0 ( r) = C 3 (r,7) , 

T 0 ( r ) = T ( r ,z) , (4.3) 

^ A 

where v 3 and T are found from the computations in the previous section. 

Substituting (4.1) into (2.1) - (2.4) and taking the 0(5) terms of 
leading order in Re yields the following linear system 



TTo + To 1 + in 


n 

r To a 


t i » 0 p . = 0 . 


(4.4) 

■ l£ + ' Sull + ;. Su = -Jp , 

To To To yM 2 


(4.5) 

* ic 7, + 1 ta = • ini 

To 2 To yM 2 „ ar ’ 


(4.6) 

- i Ct 2 c vj i fit 2 W n V] p r 

To To ” yM 2 ” 


(4.7) 

^ [- ic f + 1W 0 T + V! To r ] + fill 

(icp - 

iW 0 p) = 0 . 

L J l . 


(4.8) 


The perturbation equation of state is 


p = To p + ^ . (4.9) 

After some algebra, this system may be reduced to the following two 


equat ions 


~ 1 ~ Wo r ~ 

V,r + r Vl X% 7 c) V] 


_ 'P 

" 7M^(W 0 -c) 


To 


[ 1 + J 2 ~Z ] ' m 2 -< V ‘> 2 


(4.10) 


i& 2 ~ 1 

— (W 0 -c) Vl = - — j- p r , 
To yM » 


(4.11) 
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a rcsuli also to be found in ihc work on the siabilily of jets by Michalkc 
(1971). Writing 

Vj = Cq>. a = a/C . (4.12) 

gives 

C W on <£» 

9n + — *= — 9 - - mJL 
l+X^2 + ^ W 0 -c 


.) P 


7 M^(W o -c) 


T 0 [l + J n ^ r f H ] - M«2(W 0 -c) 2 

L a 2 (i+xc 2 + Cn) 2 J 


together wi th 

ia 2 (W 0 -c) — = * 2u_ . 

0 To yM 2 * 

These equations may be combined to eliminate p to give 


(4.13) 


(4.14) 


d 

dti 


(W °' C) [ ^ + i T Jfe ' Y ' ft 9 ^ - w ot, 9 


To 1 + 


n 2 C 2 

a 2 (1+X C 2 + £q) 2 


] • <v 


C) i 


= T5 <V c > <P • 

Alternatively 9 may be eliminated to give 

(VC)2 m2 - • T0 [ 1 * ? o r w } } 


(4.15) 


w 0T1 - } 2n 

011 1 + X& + Cq J a 2 


Pn t q_ 


(Wq-C) 


- (W 0 -c) 2_ 
dq 


T ; PH - 

a 2 (W 0 -c) 


(4.16) 


These equations reduce naturally to the system considered in I by 
setting n = 0 (i.e. axisymmetric disturbances) and X = 0 (zero cone 
angle) . 

The boundary conditions to be applied to the above systems are that 
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9 = ^ = 0 on 11 = 0, 
logcihcr wi th 

♦* » *- { K „*l ™ <«} • 

- T 9=„ M 2 „ i ay (1-c) K n (rj) 

P IT - M_2'(i'.c)2]l 

where h = ± = [1 (1 - c)2] ^ + T)J ; (4.20) 

here ihe appropriate sign is chosen to ensure that the real part of ri 
is positive as tj — * », to ensure boundedness. Equations (4.15) together 
with (4.17) and (4.18) (or (4.16) together with (4.17) and (4.19), or 
(4.13), (4.14) together with (4.17), (4.18), (4.19)) constitute an 
eigenvalue problem; here, we consider just the temporal problem, i.e. given 
n and a (real), to find c (generally complex). The problem was treated 
numerically using a Runge-Kutta method based on" (4.13) and (4.14), shooting 
inwards from some suitably large value of q, with contour indentation in 
the manner described in I. 

The ! tripiv generalised inflexion condition" 

In I, a condition for the existence of so called inviscid, neutral, 
axisymmetric mode, i.e. those for which 

1-1/M* < c < 1 + 1/M m , (4.21) 

implying disturbances that decay as q — » «, was derived. Here, we go on 
to derive the corresponding condition for non- axisymmetric modes. 

If we multiply (4.15) by <p* (where an asterisk denotes a complex 
conjugate), and subtract from the resulting equation its complex conjugate, 
then the following equation is obtained 


(4.17) 

(4.18) 

(4.19) 
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9 * d 
W 0 -c 3q 


(W 0 -c) (9q + 7 <P) * W 0T1 9 
X 


9 d f (W 0 -c*) (*,* + I 9 *) - W 0T1 9 * 

“ V* 7 

where we have written 


X = To 




M 2 „ (W 0 -c ) 2 , 


and r = (1 + X £ 2 + £q)/£. 

After some algebra, (4.22) can be written as 


d 

’ 9n + 7 <P 

d 

dq 

X 

* dq 



(4.22) 


(4.23) 

(4.24) 


= r 99 * 


1 d 

W^c aff 



W 0 -c’ 


f r v] - 

drt L x J J 


Writing c = c r + i cj , 

then the neutral state corresponds to Cj — * 0 . 

In this limit, (4.25) may be written 


(4.25) 

(4.26) 


d r [ 9* (<fci + 7 <p) * 9 ( V + 7 9* )] 

3q _ 

= 2 i r 2 1 9 1 2 c j d 
1 W 0 - c i 2 dti 

Following the arguments used in I, if the neutral mode is subsonic, 
and satisfies the impermeability boundary condition on q = 0 , then to 
avoid any contradiction the right-hand-side must be zero at the critical 
point where W 0 = c namely q = ti j ; this requires (recasting the equation 
in terms of the original variables) 




(4.27) 
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d 

3n 


Woii 


To n+x £ 2 + 



n 2 C 2 1 

a 2 ( l+X£ 2 +£q) 2 J 


= 0 . 

1) ■ Tlj 


(4.28) 


This result represents a further generalisation of the so-called 
"doubly generalised inflexion condition" (or a "triply generali-sed inflexion 
condition"!); setting n = X = 0 in (4.28) retrieves the corresponding 
result found in I. Note that in the general result for any shape of cone 
surface described by 

r = 1 + Xf(C) , (4.29) 

merely requires the "X£ 2 ” in (4.28) to be replaced by "X f(£)". 

I 

5. Numerical Results „ 

There are clearly many choices of parameters to be made in this study. 
The strategy here will be to carry out a detailed study for one choice of 
Prandtl number (0.72), ratio of specific heats (1.4), .cone angle (X = 1) 
and Mach number (M*, = 3.8) for insulated wall conditions; results are 
shown in Figs. 4-11. Results for Cj are presented for a range of values 
of n, at fixed values of £. Specifically results for £ = 0.01 (Fig. 4), 

£ = 0.05 (Fig. 5), £ = 0.1 (Fig. 6). £ = 0.2 (Fig. 7), £ = 1 (Fig. 8), 

£ = 5 (Fig. 9), £ = 20 (Fig. 10) and £ = 75 (Fig. 11). In these figures 

we just concentrate on growing/neutral modes. All results may be regarded 
as being independent of numerical grid. In a number of cases the 
distributions are shown on two figures - this is to increase the resolution 
of certain features of the distribution in regions of small growth rates. 

In all these figures, the n = 0 results are delinerated by a solidus, 
the n = 1 results by - - - -, n = 2 results by — - — - — , 

n = 3 results by — -- — -- — and n = 4 results by 

Neutral points are also marked on the axes. 
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The corresponding planar results arc shown in Mack (1965, 1987a) and 
1; the important feature in this case is the existence of two primary 
instability modes (and others, but of considerably smaller growth rate), 
with the largest growth rate being associated with the so called "second 
mode” . 

Figure 4 shows distributions at the location £ = 0.01, for 
n = 0,1, 2, 3, 4. The axisymmetric (n = 0) results are very similar to 
Mack's planar results described above (and the axisymmetric results for 
the circular cylinder case close to the leading edge, as considered in I). 
The first mode extends from close to a = 0 (where c r is approximately 
1-1/M*, up to a = 0.14 (terminating at a doubly generalised inflexion 
point); we shall refer to modes of this general type as mode _ I . A second 
mode originates at a = 0.24 (where c r = 1, and hence may be regarded as 
a generalised inflexional neutral mode, with the critical layer occuring in 
the freestream) and terminates at around a = 0.4 (which corresponds to 
a second doubly generalised inflexional neutral mode); we refer to modes 
of this general type as mode II. 

In some ways the corresponding n = 1 results for cj (also shown in 
Fig. 4) are very similar to the axisymmetric case. However a third mode is 
seen to develop, not present in the corresponding axisymmetric (and indeed 
planar) results. One important di st inct ion between the n = 0 results and 
those for n * 0 emerges in the limit as a — * 0, for which cj 4-* 0 if 
n * 0 (the limit as a — > 0, £ = 0(1) is considered in Appendix A), 

although of course the temporal growth/decay rate ctcj is nonetheless zero 
at a = 0. We shall refer to this additional mode as mode 1^. As a 
increases, mode 1^ rapidly disappears, terminating at a supersonic neutral 
point (i.e. where c r < 1-1/M,,). As a increases, a further neutral mode 
soon emerges, at a second (supersonic) neutral point; this mode is of the 
class I type described previously. Thereafter, as a is increased 
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further, the n = 1 distribution closely resembles the n = 0 results. 

The n = 2 distribution is qualitatively similar to those of n = 1; 
however when n = 3, modes 1^ and I arc seen to amalgamate, although 
mode II remains quantitatively similar to the mode II results of the 
previous n values. In the ease of n = 4, mode I admits further 
enhanced growth rates although it is still mode II that possesses the 
largest growth rate (a Cj). 

, i 

It is not surprising that at (, = 0.01 results very similar to 
corresponding planar results are obtained, since curvature, generally, will 
play a minor role in the physics here; indeed a crude examination of (4.15) 
suggests that as £ — » 0, the corresponding planar Rayleigh equation is 
attained. Howevcr.as a — » 0 and £ — i > 0 a non-uniformity is present; 
this aspect is taken up in some detail in the following section, where 
further light is shed on the additional mode 1^. 

Figure 5 shows results for c, at £ = 0.05. The axisymmetric mode 
exhibits the same general features as the corresponding mode shown on 
Fig. 4, although both modes I and II have diminished growth rates. In 
the case of n = 1. it is to be noted that mode 1^ has (already) 
amalgamated with mode I, whilst mode II has significantly reduced 

i 

growth rates compared to the £ = 0.01 value, although this remains the 
more dangerous mode. In the case of n = 2, the combined modes I and 1^ 
exhibit'an enhanced growth rate when compared with 'the £ = 0.01 station, 
whilst comparing corresponding mode IIs, we note that the growth rate is 
diminished; indeed, for n = 2, the maximum growth rate of the two modes is 
comparable. In the case of n = 3,‘ the first mode again is more dangerous 
compared to the corresponding mode at £ = 0.01, and this trend is repeated 
by the n = 4 results, although this value of n indicates that an 
increase in n is causing a decrease in growth rate. 
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Figure 6 shows results at the £ — 0.10 location. These confirm the 
trends observed previously, of a less unstable mode II, whilst the 
(combined) mode I has a maximum growth rate comparable to mode II; 
although it appears that the n = 0 mode II remains the most unstable. 

Figure 7a shows Cj versus a distribution at £ = 0.20. In this case 

mode II is barely visible, and so Fig. 7b shows these modes on a larger 
scale. Interestingly, it seems to be the n = 2 mode I that exhibits the 
largest growth rate and is consequently the most important. 

Figure 8a shows Cj distributions at £ - 1. Again, because of the very 
small growth rates, mode II is shown on an enhanced scale in Fig. 8b. 

Note that at this axial location the n = 0 and n = 4 (and above) cases 
do not possess an unstable mode I at all (whilst the n = 4 mode II has 

such small growth rates that it is not visible even on the scale of 

Fig. 8b). Overall, it is the n = 1 mode I which is the most dangerous. 

Moving further downstream to £ = 5 (Figs. 9a, 9b) we see a recovery 
in the maximum growth rate of mode II. Indeed, these results show some 
resemblance to the £ = 0.2 shown in Figs. 7a. 7b. This is not too 
surprising, given that, on account of the Mangier transformation 
(Mangier 1948, Stewartson 1964), results as £ — » - mirror those 
as £ — » 0 (except for a multiplicative factor of n/T~ in a). 

This trend is confined in Fig. 10 for £ = 20, which may be 
compared directly to the £ = 0.05 results of Fig. 5. Notice, in 
particular, the re-emergence of mode 1^ for n = 1. 

Figure 11 shows results at the furthest downstream location studied, 
namely £ = 75. The n =0 modes now correspond closely with the planar 
results of Mack (1984, 1987a for example), (with the factor VT" 
multiplying a), whilst mode 1^ is visible for n = 1 and 2, and the 
union of modes 1^ and I is clearly seen in the case of n ■ 3. 
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Thus to summarise, these results, which will guide us in certain 
asymptotic aspects in the following section, we observe the following 
general features of the stability of the flow; (i) results as £ — * ~ 
mirror those as C — * 0, except for a multiplication factor of \/3T 

(ii) Lateral curvature has a strongly stabilising influence on mode II, 

/ • . ..... 

together with the axisymmetric mode I (this is in accord with the results 
found in I); (iii) there emerges a third mode. I A , as a — » 0, with 
C — * 0 or £ — » ». in the case of n x 0. 

In the following section we go on to consider various asymptotic limits 
of the system (4.15) (or (4.16)), guided partly by the observations made 
above . 


In this section we consider a number of asymptotic limits of the 
stability problem, to give us a better understanding of the details of the 
numerical results described in the previous section. 

Perhaps the most intriging feature of these numerical results is the 

emergence of an additional mode as £ — * 0 (or £ — » ») with a ♦ 0. We 

investigate this feature first. 

— C -=± Q. — g-= flip q i C — » q = Q (C~h 

Since the problem as posed is basically equivalent as (j — » 0 and 
C — » we consider just the former limit, and later we describe, briefly, 

how the results for the latter limit can be simply inferred. 

As noted in Section 5, as £ 0, (4.15) is seen to generally reduce 

to the planar system treated by Mack (1984, 1987a for example). However this 
will no longer be the case if a = 0(£), since then the denominator on the 
left -hand- side of (4.15) no longer reduces to the planar result. 
Specifically, let us write (consistent with (4.12)) 
a = Ca , 


( 6 . 1 ) 
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where ii is assumed 4 = 0(1) as £ — » 0. The results for mode 1^ 

shown in the previous section, together with other results obtained by the 

authors indicate that as £ — ♦ 0, c — » 0 also. 

Partly guided by this, if q = 0(1) we set 

c = £ c i + ^ 2 c 2 + C^ c 3 + • (6.2) 

9 = 9o( T i) + C9lOl) + C 2 92< 1 1) + C 3 93(^) + ( 6 - 3 ) 

W o = W 00 (ti) + C W ol Oi) + ; 2 W o2 (ti) + C 3 W o3 (0 + •••■ (6-4) 

T 0 = T 00 (q) + £ T 0 j(q) + C 2 T'o2( 1 l) + C 3 ' r o3( T l) + ••• . 

(6.5) 

where W 00 (tj ) and T 00 (r| ) represent the planar values of the velocity 
and temperature profiles repectively, and WQi(q) and Toi(q) etc. 
correspond to the perturbations to the basic flow caused through curvature. 
To leading order, equation (4.15) reduces to 


d 


w oo 9oq • 

T 0O [1 + n 2 ] 

ft* 


W OOTl 9o 


- w: 


00 



( 6 . 6 ) 


i . e . W 00 9ori ’ ^ooq 9o “ ^o 



(6.7) 


where k 0 is a constant. However we require that 9o( T l=0) = 0* whilst 
<p 0 is restricted not to grow exponentially as q — ♦ Consequently we 

must have that 

k 0 = 0 . (6.8) 

and hence 

9o = A o w oo0l) • < 6 - 9 > 

where A 0 is an arbitrary constant (i.e. the unknown amplitude of the 

eigensolut ion) . 

Equation (6.6) is not a uniformly valid approximation to (4.15), for all 
q; specifically a breakdown occurs when q = 0(£ _1 ). We define 


q = 1 + £q = 0(1) , 


( 6 . 10 ) 
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(i.c. q represents a scale comparable to the radius of the cone), and 
on this scale we expand 9 as follows 

<9 = * 0 (n) + C (n) + ... • (6.11) 

To leading order equation (4.15) reduces to 


d 



1 



n 2 1 

521? J 


- M ,,, 2 



( 6 . 12 ) 


The solution to this, which matches on to (6.9) as f| 


*o - 


K n [i a (Mo . 2 - l)jjn 

K„* [i a (M » 2 - 1)*] 


i s 

(6.13) 


where K n ( z j ) is the Bessel function of order n, argument zj (the 
K n (zj) solution is chosen in preference to the Im (zj) solution in 
order that disturbances are propagated along characteristics in the 
downstream direction - see Ward 1955, Kluwick et al 1984, Duck and Hall 
1989a, b). 

Returning now to the q = 0(1) layer, the 0(£) correction to 9 is 
given by 


d_ 

dq 


Wpl ^pq - ci 9 0 g + W 00 tpiTj + W 00 9 0 . Wqqh 9j - Wq]^ f 0 
T 00 11 + n 2 ] - M ~ 2 w oo 2 

5? 

(6.14) 


and so integrating this equation once we obtain 


Wol 9oq ' C 1 9oq + w oo 9lq + w oo 9o * ^ooq 9l * w olq 9o 


= kl 


l 00 


[1 + 


n 2 ] - M .. 2 W c - 2 


"00 


(6.15) 


where kj is an arbitrary constant. 9 j must not be exponentially 
large as q — » «, and consequently we must have 
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in 

together with 


+ A 0 - k. 


1 - M 2 + n 2 


(6.16) 


- cj A 0 W 00 q(i, - 0) - kj T oo (0) 

However in order to match correctly with (6.13) 


1 + n 2 

5* 


(6.17) 


<PlT, 


= ^ U = 1 (6.18) 

= A n itt (jQ . p* K n ' ' [i& (M „ 2 . \ji] 

K n ' lie t (M.2 - l)il 


Eliminating kj , we obtain the following result for c j 


Cj = 


l + i a (MqJ - i ) 1 K n * * [ia qaJ - i)j_] 

K n '.[ift (M « 2 - D*] 


(6.19) 


Too (T1 _ 0) 


1 + n 2 

a 2 


Woor, 01 = 0 ) mJ - 1 - n 2 ] 

a? 

The asymptotic forms for this expression in the limit of large and small 
a may be found readily. Firstly as a — > » we have 


C! -r — T QQ (T1 ~ 0) r + Ofa' 1 ) 

Wool, (T1=0)(M«2 . i)l 

Secondly as a — * 0 

c ' - "w!s:; iSi + °< a2 > for n *' • 


(6.20) 


T oor, 


c i - • w!°S ' j^g T + ° <a2 ioga) for 


( 6 . 21 a) 

( 6 . 21 b) 


’OOT, 

In fact it is quite easy to show that as a — » 0. 

j , , - a 2n (M- 2 - l) n Too ( 0 ) 2 2 2n 

Im ci} — ♦ = >5 . ( 6 . 21 c) 

1 Woot, ( 0 ) . [ (n- 1 ) ! ] 2 

Equations (6.21) are precisely the (real) values found by Duck and Hall 
(1989b) for the downstream limit of a non-axi symmet ric viscous mode 
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(taking into account the different scalings used in Duck and Hall's paper). 
Consequently we expect that as £ — » 0/«. on a scale smal ler/1 arger than 
that of the cone radius, we expect this mode to become predominantly viscous 
in nature, and to be described by triple-deck theory. 

Distributions of I m { c j } with ft are shown in Fig. 12 for n = 1,2,3 
(with (6.20) also shown). Unfortunately (perhaps) it is seen that 
Im(ci) < 0 for all a (confirmed by (6.20)). From (6.21) we also have 
that Im{ c j ) = 0(0 as £ — » 0. Unfortunately, also, the 0({|2) and higher 
correction to this mode would require a large amount of algebra. However we 
are able to make progress, in particular obtain an estimate for 
Im(c (a = 0)) by considering instead the limit as £ — » 0 of equation 

(A. 2), pertainent to the a — * 0 case. (With a = o(0 *we consider this 

aspect in the following subsection). 

Finally, for this subsection, note that the £ — * ~ results may be 
simply inferred from these £ — » 0 results, simply be replacing the small 
parameter in the various expansions by the small parameter '1A£'. 

More subtle differences between the £ < < 1 and £ > > 1 solutions only 
appear at higher orders. 

fL 2 fl LsJ L C — t Q (fil ( 

The system (A. 2) turns out to be rather easier to analyse as £ — * 0 
than does the corresponding finite tt (= ct£’S problem. We again utilise 
expansions (6.2) - (6.5) (although see (6.48) for the q> expansion). 

To leading order, we have for q = 0(1) that 

q >0 = A o Woo(n) < (6.22) 

where A 0 is some (arbitrary) amplitude parameter. 

At the next order we have the following system 

^oo ‘Plq ' A o ^oor| C 1 + A o ^oo^ 

*^oor| ^1 + A o W 0 1 ^ooq * A o ^olt| ^oo. = ^1 ^oo • (6.23) 
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where kj is a constant, and we have utilised (6.22). Selling T| = 0 in 
(6.23), assuming (^-q = 0 , then 

- A o Cl W OOT ,01=0) = k! T 00 (q=0) . (6.24) 

The boundary conditions as q — > «• must be compatible with (A. 3), together 
with (6.22). Defining 

t = 1 + n/C = 0 ( 1 ) . (6.25) 

then we must have an outer solution of the form 

<P° Ut = A r'"' 1 . (6.26a) 


where A = A 0 + £A] + (6.26b) 

and so we must also have 

I out | 


9m 


- 9or 


r = j = - (n+1 )A 0 


(6.27) 


Substituting this into (6.23) yields 


ci = 


_ n T*no (jj ~ 0) 


(6.28) 


*oot| m - o) 

in accord with (6.21). In order to estimate complex values of c we must 
proceed to higher orders in £. 

At the next order in £ we obtain the following equation governing 92 


w oo 92t| + 2X W 00 9 0T i + 2q W 00 <p lTJ + W o2 9 0T) 

*^o2q 9o + ^ol 9lri + ^ w oo 9 oti ' C 1 9lq 

t 

+2W 0l n <p 0T , + W ol 9 0 - W oln 91 * 2rj W olll 9 0 
■2 tj c j 9 0T | - C 2 9otj + ^00 9l + *1 ^oo 9o 
* C 1 9o " ^ooti 92 ’ 2X W 00T j 9 0 ■ 2t| 9 j W 00T j 

^oot] 9 o = ^2 ^00 + ^00 | — dT| 

0 

+ ki T ol . (6.29) 


We shall defer any consideration of this equation, and move to the next 
order of which yields 
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^oo 93 t| + 2 ^oo 91 ti + 2 ^l^oo 92 t| + 2 Xtj Wq 0 9oq 

+ n 2 W 00 (p ln - C! <p 2l1 - 2Xc) (p 0T1 - 2i] cj 9lq 

- n 2 Cl <Poti - c 2 <PlTi - 2r \ c 2 9011 - c 3 9oii 
+ W 00 92 + ^ *oo 9o + ^ W 00 9l * C 1 91 

* Cj T1 <p 0 - c 2 9o ■ W 00T1 93 • 2 X w oon 91 

- *1 2 W ohl <p 0 - 2 t| 92 W 00T1 - q 2 W 00T j 9j 

+ W ol 9 2 ^ + 2X W ol 9 0T1 + 2r) W 0 i 9j^ 

+ ti 2 W 0 i 9 0T j + W o2 9]^ + 2 ti W o2 9 0T j 
+ W 0 3 9 0T | + W 0 j 9 j + W o] T] 9 0 + W o2 9o 
'* w o3ti 9o * 2 n W 0 2 tj 9o " 2 ^ ^olT| 9 o 

- 2Xt) W 00T | 9 0 - W o2 -q 9j - 2 t] W 0 ]t| 91 

* w o1ti 92 = k 2 T 0 i + k i T o2 + k3 T 00 


n 2 T, 


f9l T oo ' 9o T 0 i] W 00 ' C 1 Too 9o + *ol 9o T 0 o 


oo 


dii 


>00 


+ n 2 T, 


ol 


9o W 00 


(6.30) 


‘00 


Our main goal here is now to determine the leading order imagiDiiy 
component of the complex wavespeed c (we do, of course, already know the 
leading order real term). Now since the above equations just contain real 
coefficients, any imaginaries must, of necessity, only arise at a critical 
point, where, c = W 00 . Since c = 0(0. this must occur when H = 0(C). 
We therefore consider a thin layer relative to the r\ = 0(1) scale, namely 


tl = H/C = 0(1) . (6.31) 

On this scale, the expansion for 9 is expected to develop as 

9 = Oo 01> + C 2 9j (tl) + C 3 9 2 (tl) + (6.32) 

where the Oj are expected to be normal i sed in such a way as to be 
generally 0 ( 1 ) quantities. 
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It is easy to show that 


4»o(tl) = Aq W 00l] (ti= 0) t) 

(where A 0 was introduced in (6.22)), and also 


(6.33) 


*i 01) = Aj W 00T1 (11=0) h . (6.34) 

where A j is a constant, linearly related to A 0 , and we have used the 
property that W 00T ^ (ti=0) = 0). 

After some algebra, the equation for ^(’l) ma y be wr » ttcn 
[W 0 or, (0) ^ - c,] 4> 2 q * w oon (°> *2 

= k 2 + 3i) ci A 0 W 00T) (0) + ^ r\ A 0 (0) W OOt> (0) 

- [Woor, (0)] 2 A 0 r\ 2 

-2 

. w oon (onwffi C 1 A 0 JL , (6.35) 

where here k 2 is an (arbitrary) constant. 

If we take (as we are quite at liberty to do) A 0 and k 2 to be l&a± 
constants (this is not essential for our arguments, but simplifies the 
following argument), then we now consider just 4>2* (where here and 
elsewhere a superscript i denotes an imaginary component). This quantity 
is triggered by the well known + i re jump in the logarithm (Mack, 1984 for 
example) across the critical layer. Specifically, here, this is caused by 
the q dependancy on the right -hand-side of equation (6.35) (k 2 plays no 

role in this). If (6.35) is written symbolically as 


[Woon (0) - ci] <I>2^ - W 00T1 (0) d>2 = R 


(6.36) 


then 


*2 = [ w ooq (0) V * ci] 


R dq 


71 • 


(6.37) 


[Wooti (Q) ^ * c l^ 

Evaluating this integral, taking just the imaginaries together with the 


limit as tj 


00 


yields 
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*2 ~ K A o [ w oo n (°) ti c,] cj 


+ c > 


W ooqqq<°> 

[ w OOll ( 0)] 3 


‘00 


Tin 


(0) Cj 


T OO (0) [^ooq(O)] 


= b' [w oon (0) ti - cj] . 


(6.38) 


where B 1 = 7tA 0 Cj 


+ Cj 


’00 


nnn 


(0) 


Wqoii(O) [w 00 ll co)]- 


r 00TlT1 (0) C 1 


(6.39) 


^ 00 ( 0 ) [Woo^W] 

Equation (6.38) then provides the lower boundary conditions for the system 
(6.29) and (6.30). 

Since (6.29) contains just real coefficients (and taking C 2 to be 
real, an assumption that may be justified a posteriori), then we must have 

<P2 = B‘ W 00 (r|) , (6.40) 

where B 1 is defined in (6.39). 

If we take the imaginary part of (6.30) and allow tj — * «, then 

* 3 n In -- - *2 I , - - - ^ • < 6 - 41 > 

However we require (on account of (A. 3)) that 
i 


9 3qlq - - ~ - (n+1) 9 2 1 |ti 


(n+1) B 


k 3 ‘ = - nB 1 


and so 

Setting q = 0 in (6.30) yields 


(6.41) 

(6.42) 


‘ C 1 9 2q |q = 0 ’ c 3 A ° W <W°) * w ooq(°) ^3 |q = 0 

i (6.43) 

- - nB Tq O (0) 
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However from (6.38) we have 

* 21 / |q = 0 = B ' W OOT1 (0) • 


(6.44) 


9 


3 h = 0 


= c j B' 


(6.45) 


Consequently, (after substituting for Cj), 


. T (0) Jt 

1 2 00 

c = n j 

-J r«i / r \ \ T *» 


[Wooq(O)] 


" T oo(°> w oowi (0) 


^OOqW 


[Wooq] 


W^ n l 

[ w ooq(°)] 

(6.46) 


Using the governing equations for the basic flow allows this expression to 
be simplified slightly, to 


" 2 T oo <°> 2 « 

[*001] W] 


1 + 2o M 2(y- 1 )n 

OP 


(6.47) 


In fact the expansion for <p in (6.5) is not quite complete as it 
stands, since the analysis of the 1) = 0(1) layer above indicates the 
presence of loga-rithmic terms; specifically we require 
9 = 90^) + ?9i0l) + C 2 92( T l) + C 3 93( | l) +• • • 

+ log C [ C 2 921 01) + C 3 931 Ol> ] • (6-48) 

where 921^) = A 21 Woo^) - (6.49) 

wi th A 21 a constant . 

A comparison of the fully numerical computation of Real j c (a — 0) J 
with the asymptotic formula (6.28), as £ — » 0 is shown in Fig. 13. The 
agreement is seen to be entirely satisfactory. Unfortunately the 
correlation between the computed Im j c (a = 0) J and that obtained 
using (6.47) is much less satisfactory. However, this poor correlation 
is not unexpected for two reasons. Firstly accurate computations of 
Im | c J in this limit become exceedingly difficult, as confirmed by the 
quite complex asymptotic structure detailed above, with both thin 
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(il = 0 (£)) and thick (rj = 0 (I/O) lcngthscalcs emerging. 

Secondly the asymptotic form for Im j cj J is achieved very slowly as 
C, — * 0, at least in one particular configuration, where, with n = 1, the 
imaginary wavespeed has a leading order coefficient of approximately 
3.898 x 105 £3, A comparison between the numerical and asymptotic results 
is not shown in this case. 

In the case of £—♦«>, wc ma y j usl replace the small parameter 'O 
in the above, by the small parameter '1MO- In the following subsection, 
we consider the behaviour of mode I, as ( -* 0. 

6.3 i -t d ct = .Q(Ctl 

The numerical results presented in the previous section strongly 
suggest that for the most part, as £ — ♦ 0 mode I has features very 
similar to the planar case, for all values of n. However there is one 
important exception to this, namely the behaviour of the lower neutral point 
in this limit. The planar case has that c — * 1 - 1/M., as a — » 0, 
corresponding to the so-called "sonic" mode. However in the case of our 
numerical results, there is evidence of a shift in this neutral point, 
along the positive real axis, and the neutral point becomes (slightly) 
supersonic, with c < 1 - 1/M«, as £ — » 0. 

A (sensible) balancing of terms suggests that we might look for a 
solution of the form 


wi th 


C — Cq + c j + . . . 

9 » <Po + C<Pl +• •■ 

W 0 = W 00 + c W 0 1 +.. 

T 0 = T 00 + (, T 0 ] +. . 


a = a, a = 0( 1 ) 


(6.50) 


(6.51) 
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To leading order, (4,15) yields 


d 

3n 


(W 


00 


c ) <p - W <p 

o' y OTJ 0T| t 0 


= 0 


(6.52) 


where 

*0 = T oo - M ~ 2 ( w 00 - c 0 ) 2 . ( 6 - 53 ) 

The solution to (6.52) which satisfies the impermeability condition on 

t| = 0 is 


9o = K 0 < W oo ’ c o) T 


,dn . 


(6.53) 


(Wqo-Co) 4 - 

V 

where K 0 is some arbitrary constant, and the integral is to be taken 
underneath the critical point. 

Since we require that <p 0 — » constant as q — » •* (in order to match 
with (4.18)) we must have that 

c 0 = 1 ± 1/VL . (6.54) 

in order that the integral (6.53) remains bounded as q — » Further, we 

take the ne ga t i ve sign to be consistent with the numerical results and 
or comments above; indeed, this is simply a repeat of the planar calculation 
(Lees and Lin 1946) . 

Curvature plays an important role at next order, namely 0(£). The 
governing equation in this case is 


d 

3q 


^ ^ (Woo ' c o)( t PlT| + Vo) " C 1 9oq ’ W 0 oq 9l + Wol 9oq * Wolq 9o j 

£ ( w oo ’ c o) 9oq ‘ W 0 oq Vo j ^ + 

+ 2M„2 (W 00 - c 0 ) C! - (W 00 - c 0 )W 0 i jj 

_ ft 2 (Wqo ~ c q) .9q _ (6.55) 

i nrv 


Consequently, using (6.53), 
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(Woo " c 0 ) (<Plri + 9o) ' C 1 9oii ‘ W 0 or| 9l * W 0 i-q <Po 

X 

K 0 [ + T ol + 2M J (W 00 - c 0 ) Cl - 2K? (W co - c 0 ) W ol J 


- Kj T 0 + T 0 


(W - c ) <p 
v 00 0 Y o 


•oo 


dii , 


(6.56) 


0 


where Kj i s a constant . 

In order to match correctly at the outer edge (q — * ■») with (4.18), 
we must have that 


9 


IT] 


% K n <V K o 1 
K n ' (H 0 ) M- 


(6.57) 


where 


i = + 


dri , 


(6.58) 


6 (Wqo * c 0>' 


and T)o = mJ a (2c 1 ) i . (6.59) 

Talcing the 0(1) terms of (6.56) as q — » •» and using (6.57) yields 
the following non-linear dispersion relationship for c j : 


Ho ^n Olo> , j . I" n2 + ou Cl 1 

C'T & + 1 - -r [ p + 2M “ c i J • 


(6.60) 


The integral (6.58) was evaluated numerically, and for the conditions 
prevailing in all our numerical results it was found I - -228. 4-59. 3i . 
Equation (6.60) was solved using Newton iteration, and results for 
Real ( c i } an d Im j cj J fof various n are shown in Figs. 13a and 13b 
respectively. 

Notice that as a — » <», (6.60) predicts that one family of solutions 
has the property that 

ci - ; 2 i 2 
2 MS 


(6.61) 
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which is in agreement with the a « 1 expansion for c for the planar 
ease, namely 

c = 1-1/M„ + + 0(a4) (6.62) 

2 

Equation (6.61) is shown as an asymptotic on Figs. 14a and 14b. Note 
that the (real) family of cj which also exist may be found as an exact 
solution to (6.60), namely 

C 1 = - , (6.63) 

2 «2 

although the importance of this mode is not thought to be great. The 
complex families of cj's are seen to terminate at a finite value of a, 
corresponding to the (lower) neutral point of mode I. Notice that in all 
cases because Real { cj J < 0 at the termination point, these modes 

I 

correspond to supersonic neutral modes. 

From the result shown in Fig. 14b, we are therefore able to offer an 
estimate for the position of the lower neutral point of mode I as £ — » 0. 
In particular, for the freestream conditions considered throughout this 
paper, for n = 0 this position is given by a.~ 0.1 0 , for n = 1 by 
a = 0.20 £l, and for n = 2 by a = 0.295 0 . Comparing these asymptotic 
results with figure 4 in particular reveals a fair degree of agreement. 

In the case of C, » 1, the above results may be easily transposed, 
by the replacement of 'O by '1A£'; the corresponding positions for 
lower neutral point are then a - 0.1 AO * for n = 0, 
a - 0.20 AO for n = 1, and a ~ 0.295 AO * for n = 2. These 
results are seen to agree quite well with the £ = 75 results shown in 
Fig. 11. 
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1. Conclusions 

In ihis paper we have studied the supersonic boundary layer flow, and 
the inviscid stability thereof, over a sharp cone with adiabatic wall 
conditions. The basic flow is seen to evolve from one planar state, to a 
second, as predicted by the Mangier (1946) Transform. 

The "triply generalised inflexion condition” is derived, this being the 
necessary condition for subsonic neutral modes, and is a (second) 
generalisation of the well known generalised inflexion condition. 

Significantly, the numerical results point to the occurence of a third 
mode of instability, not found in similar planar and (more recently) 
axisymmctric studies. An asymptotic study of this mode shows this mode to 
be linked to a viscous mode found by Duck and Hall (1989b), a study based on 
triple-deck- theory. 

The "sonic” neutral mode is found to be altered by curvature (and in 
fact becomes a supersonic neutral mode as revealed by the asymptotic 
analysis valid as C — » 0 and £ — » ».) 

Significantly, our results show that the so called second mode (mode II) 
is not always the most unstable, at least in the case of non-axi symmet ric 


modes . 
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Appendix A The zero number limit 

In the case of axisymmctric modes (and indeed of planar modes, as 
considered by Lees and Lin 1946, Mack 1965a, b, 1984, 1987 for example), 
as a — » 0 . the wavespecd c approaches the sonic limit, i.e. 

c — > 1 ± 1/M„. (A . 1 ) 

In the case of non-axi symme t r ic modes, however, this is no longer the case. 
As clearly illustrated in the results, when n * 0, Cj •/-> 0 as a — » 0. 
The explanation is as follows. 

If we (simply) allow a — » 0 (assuming n * 0), then (4.15) reduces 
to 


d f (w 0 -c) [1+AC 2 + Cn] [(l+xc 2 + Cn) <&, + C<p] 

dt| T 0 


- ' v on[ 1+x C 2 + ^] 2 <p 
T o 

= - n 4 lj (W o- c > (A. 2) 

4 0 

n 

As q — » «, this system clearly admits solutions of the form 

9 ~ T|~ n " 1 (A. 3) 

(together with 9=0 on q = 0 ) , which is completely compatible with 

. A A 

the outer solution, where q = 0(1), (ti defined by (4.19)) namely (4.18). 

Equation (A. 2) then represents a reduced problem as „ a — » 0, and indeed 
a reduced form of the triply generalised inflexion point condition (4.28) 
exists in this limit, namely 


d [ W 071 ( 1 +AC 2 + Cq ) 2 
3q • T 0 


(A. 4) 


The system (A. 2) - (A. 3) was solved in a number of cases (in an 
identical manner to the a = 0 ( 1 ) eigenvalue system) and its correctness 
was confirmed. Notice, however, since the actual temporal growth rate is 
acj, then this still reduces to zero as a — > 0 . 
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Fig. 1 Layout 
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g.7a Distributions of cj with 
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Fi g . 7b Details as Fig. 7a, enlarge 







